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Abstract 

Given a Gibbs point process P* on R'* having a weak enough potential 'if, we consider the 
random measures /ix ■= Y.x&v'^'nQ^ ^{x,V'^ n Q\)d^/xi/d, where Qx := [-A^'''*/2, is 
the volume A cube and where £,{■,■) is a translation invariant stabilizing functional. Subject 
to '5 satisfying a localization property and translation invariance, we establish weak laws 
of large numbers for X~^fix{f), f a bounded test function on R'', and weak convergence of 
\~^^^ fixif), suitably centered, to a Gaussian field acting on bounded test functions. The result 
yields limit laws for geometric functionals on Gibbs point processes including the Strauss and 
area interaction point processes as well as more general point processes defined by the Widom- 
Rowlinson and hard-core model. We provide applications to random sequential packing on 
Gibbsian input, to functionals of Euclidean graphs, networks, and percolation models on 
Gibbsian input, and to quantization via Gibbsian input. 

1 Introduction 

Functionals of large complex geometric structures often consist of sums of spatially dependent 
terms admitting the representation 
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(1.1) 



where Af c M'' is locally finite and where the function ^, defined on all pairs {x,X), with x & X, 
represents the interaction of x with respect to X. When X is a random n point set in R'^ (i.e. a 
finite spatial point process), the asymptotic analysis of the suitably scaled sums (1.1) as n ^ oo 
can often be handled by M-dependent methods, ergodic theory, or mixing methods. However there 
are situations where these classical methods are either not directly applicable, do not give explicit 
asymptotics in terms of underlying geometry and point densities, or do not easily yield explicit 
rates of convergence. Stabilization methods originating in [23] and further developed in [3, 24, 26], 
provide another approach for handling sums of spatially dependent terms. 

There are several similar definitions of stabilization, but the essence is captured by the notion 
of stabilization of the functional ^ with respect to a rate r > homogeneous Poisson point process 
V := Vt on R'^, defined as follows. Say that ^ is translation invariant if ^{x, X) = ^{x + z,X + z) 
for all z G M**. Let Br{x) denote the Euclidean ball centered at x with radius r G M+ := [0, oo). 
Letting denote the origin of W^, we say that a translation invariant ^ is stabilizing onV = Vr if 
there exists an a.s. finite random variable R := R^{t) (a 'radius of stabilization') such that 

^(0, V n Br{o)) = ao, V n Br{x) u A) (1.2) 

for all locally finite .A C R'' \ Br{0). 
Consider the point measures 

where 6x denotes the unit Dirac point mass at x whereas Q\ := [— A^/'*/2, is the A- volume 

cube. Let B{Qi) denote the class of all bounded / : Qi ^ R and for all random point measures fj. 
on R"^ let (/, /i) := / fdn and let /z := /x — E [fi]. 

Stabilization of translation invariant ^ on V, as defined in (1.2), together with stabilization 
of ^ on n Qa, A > 1, when combined with appropriate moment conditions on ^, yields for all 
/ G B{Qi) the law of large numbers [22, 25] 

lim A-^(/,ma) = rE[^(0,P)] / f{x)dx in and in (1.4) 

and, if the stabilization radii on V and V n Q\, A > 1, decay exponentially fast, then [3, 21] 

lim A-iVar[(/, /XA)]=ry«(T) / f{x)dx, (1.5) 

where for all r > 

V^{t) := E [^(0, Pf] +T [ [E^(0, P U {z^z, V U {0}) - E [^(0, Vf\\dz. 
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Additionally, the finite-dimensional distributions ((/i, A ^/^T*^), . . . , (/fc, A ^^^/S^)), /i, • • • , /fe G 
converge to a Gaussian field with covariance kernel 



The limits (1.4)- (1.6) establish asymptotics for functionals and measures defined in terms of 
independent input and one might expect analogous asymptotics for functionals of dependent input 
subject to weak long range dependence conditions. The main purpose of this paper is to show 
that this is indeed the case. We establish the analogs of (1.4)-(1.6) when V = Vt is replaced by a 
weak Gibbsian modification having an exponentially localized potential; see Theorems 3.1-3.3 for 
a precise statement of the limit theory for functionals of Gibbsian input. Gibbsian point processes 
covered by this generalization include, for low enough reference intensity r, the Strauss process, the 
area interaction process, as well as point processes defined by the continuum Widom Rowlinson and 
hard-core models. Gibbsian point processes considered here are intrinsically algorithmic. Their 
computational efficiency yields numerical estimates for asymptotic limits appearing in our main 
results. 

Functionals of geometric graphs over Gibbsian input on large cubes, as well as functionals of 
random sequential packing models defined by Gibbsian input on large cubes, consequently satisfy 
weak laws of large numbers and central limit theorems as the cube size tends to infinity. The precise 
limit theorems are provided in sections 6 and 7, which also includes asymptotics for functionals of 
communication networks and continuum percolation models over Gibbsian point sets, as well as 
asymptotics for the distortion error arising in Gibbsian quantization of probability measures. 

2 Gibbs point processes and their stabilizing functionals 
2.1 Gibbs point processes with localized potentials 

Throughout 5" denotes a translation invariant functional defined on locally finite collections of 
points X in M'' and admitting values in R+ U {+cxd}. By translation invariant we mean ^{X) = 
^{y + X) for all y G . In the sequel we refer to ^I* as the potential, Hamiltonian or energy 
functional, with all three terms used interchangeably. For a locally finite point set X in W^' and 
an open bounded set _D C M'' we define £,{X) := '^{X n D). We shall always assume that for all 
open, boimded D the potential d{"P) admits finite values with non-zero probability, where we 
recall that V :=Vt is the Poisson point process of some arbitrary but fixed intensity r > in R''. 




(1.6) 
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Moreover, we assume the Hamiltonian is hereditary, that is to say if ^{X) = +00 for some X then 
'^{y) = +00 for all y ^ X. This puts us in a position to define the Gibbs point process given 
in law by 

dC{V) ^ ' Z[^n] ' ^ ' 

where Zl^o] '■= ]Eexp(— ^^^(■p)) is the normalizing constant for (2.1), also called the partition 
function. 

The following definition is central to this paper. 

Definition 2.1 For a decreasing right- continuous function -0 : — » [0, 1] with limr^oo V'(^) = 
we say that the Hamiltonian is ^-localized iff for each x £ R'^, each finite X and each r > the 
add-one potential, inheriting from 4* its translation invariance property, 

A{x, X) ■= ^{X U {x}) - ^{X), 

satisfies 

Q < Ay^^{x,X C\Br{x)) < A{x,X) < A^''\x,X nBr{x)), (2.2) 

where A[r](-,-) and AM(-,-) are certain translation invariant deterministic functionals such that 
uniformly m a; G M'' and X cR'^ we have 

< exp{-A[^]{x,X n Brix))) - cxp{-A^'\x, X n Br{x))) <^{r). (2.3) 

In other words, even though determining exactly the value of the add-one potential A(a:, X) 
may require the knowledge of the whole configuration X, knowing just XnBr{x) we can determine 
the value of exp(— A(a;, X)) with accuracy at least tp{r) which tends to as r — !■ cxj. In case where 
both \E'(A' U {x}) and ^'(A') are +00 we set by convention A{x,X) := 0. We also require that 
A{x, 0) < +00 to prevent the Gibbs process "P* from concentrating on 0. The functionals A[r](-, •) 
and A[''1(-, •) will be called lower and upper add-one potentials respectively. Note that the required 
non-negativity of the add-one potential is not particularly restrictive because whenever the add-one 
potential admits a finite lower bound, possibly negative —a < 0, it can be reduced to the present 
setting by adding a\X\ to 5* and by replacing the underlying intensity t with rcxp(o). Imposing 
the presence of a lower bound for the add-one potential or other related growth conditions is a 
usual assumption to avoid density explosions and infinite values of the partition function in (2.1), 
see [28]. 

Every Poisson point process has a ■(/'-localized potential, since in this case = and thus A = 0. 
Less trivially, a large number of Gibbs point processes, including those in modelling problems in 
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statistical mechanics, communication networks, and biology have '^-localized potentials. This 
includes the Strauss process, the area interaction process, processes having finite and infinite range 
pair potential functions, and the hard-core and Widom-Rowlinson models; see section 5 for details. 

2.2 Graphical construction of Gibbs point processes with locahzed po- 
tentials 

For a ^/'-localized potential ^ the resulting Gibbs point process "P* admits a particularly convenient 
graphical construction in the spirit of Fernandez, Ferrari and Garcia [10]-[12]. While adding a 
number of new ideas, in our presentation below we follow [10]-[12] as well as the developments 
in [4]. Consider a stationary homogeneous free birth and death process (p^)(gR in D with the 
following dynamics: 

• A new point x G D is born in pp during the time interval [t — dt, t] with probability rdxdt, 

• An existing point x G dies during the time interval [t — dt, t] with probability dt, that is 
the lifetimes of points of the process are independent standard exponential. 

Clearly, the unique stationary and reversible measure for this process is just the law of the Poisson 
point process V f\D. 

Consider now the following trimming procedure performed on p^ , based on the ideas developed 
in [10]- [12]. Choose a birth site for a point x G D at some time i e R and draw a random 
number vj G M+ from the law given by the distribution function 1 — '(/'(■)■ Then, accept it with 
probability exp(— AM (a;, n Br,{x)))/ilj{r]) and reject with the complementary probability if 
the acceptance/rejection statuses of all points in pP_ n B^{x) are determined, otherwise proceed 
recursively to determine the statuses of points in Bj^{x). 

Before discussing any further properties of this procedure, we have to ensure first that it actually 
terminates. To this end, note that each point x with the property of having the ball B^(x) devoid 
of points from pP_ at its birth time t has its acceptance status determined. More generally, the 
acceptance status of a point x at its birth time t only depends on the status of points in pP_r\B^{x), 
that is to say points born before and falling into Bn{x). We call these points causal ancestors of x 
and, in general, for a subset A <Z Dhy Ant [A] we denote the set of all points in n A. their causal 
ancestors, the causal ancestors of their ancestors and so forth throughout all past generations. The 
set Ant [A] is referred to as the causal ancestor cone or causal ancestor clan of A with respect to 
the birth and death process (pf )teR. 
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It is now clear that in order for our recursive status determination procedure to terminate for 
all points of pf in A, it is enough to have the causal ancestor cone Ant[A] finite. This is easily 
checked to be a.s. the case for each A C D - indeed, since D is bounded, a.s. there exists some 
s <t such that = % and thus no ancestor clan of a point alive at time t can go past s backwards 
in time. 

Having defined the trimming procedure above, we recursively remove from p^ the points re- 
jected at their birth, and we write (7/')teM for the resulting process. Clearly, 7/^ is stationary 
because so was pf and the acceptance/rejection procedure is time-invariant as well. Moreover, the 
process 7/^ is easily seen to evolve according to the following dynamics: 

• Add a new point x with intensity t exp{—A{x,^f j)dxdt, 

• Remove an existing point with intensity dt. 

These are the standard Monte-Carlo dynamics for (ID as given in (2.1) and the law of is 
its unique invariant distribution. Consequently, in full analogy with [10]- [12] the point process 7/^ 
coincides in law with for alH G M. 

2.3 Exponentially localized potentials and infinite volume limits 

Recalling the definition of ip from Definition 2.1, we henceforth assume that there is a Ci > such 
that 

V'(r) < exp(-Cir) Vr > 0. (2.4) 

It should be emphasized that we require (2.4) to hold for allr > and not just for r large enough. 
It is known, see [10]- [12] where a proof based on subcritical branching process domination is given, 
that if Ci is chosen large enough, then all causal ancestor cones are a.s. finite and, in fact, there 
is a C2 > such that for all t,Re M+ and ^ C .D we have the crucial bound 

P[diamAnt[^] > R + diam(A)] < vol(^) exp(-C2i?). (2.5) 

Moreover, the constant C2 in (2.5) does not depend on D. If (2.4) is satisfied with the constant Ci 
large enough so that (2.5) holds as well, then the potential \1/ is declared exponentially localized. 

Putting D„ ■= [— n. n]'^, this puts us in a position to construct the infinite volume limit (ther- 
modynamic limit) for as n — > 00. Indeed, consider the infinite volume version pt of our 
stationary free birth and death process , constructed as p^ with D replaced by K'^. Clearly, 
for each t €R we have that pt coincides in law with V. Moreover, in view of (2.5) and recalling 
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that c there did not depend on D, we see that the trimming procedure as described above is also 
vaHd for the infinite volume process pt, yielding the stationary trimmed process 74. These remarks 
justiiy defining the following point process, used in all that follows. 

Definition 2.2 We define the thermodynamic limit := to be the point process coinciding 
in law with 70 and hence with 7* for all t. 

To provide some further motivation for granting to the name of thermodynamic limit note 
that the process enjoys the following important property: for any bounded set D C M"^, any 
locally finite point configuration X C and any finite point configuration y C D the conditional 
law of n D on the event D D'^ = X is given by 

n n = Af) _ exp(-^(y| A-)) 

dC{VnD) ZdI^X] ' ^ ' 

where 

= Eexp(-*(P n D\X)) 

whereas 

<i>{y\x) := lim ^{y UXn B{0, r)) - n B{0, r)) 

r — *oo 

with the existence of the limit guaranteed by the localization condition (2.2). Moreover, the 
so-constructed T'* is the only point process with the above properties - to see it take Dn '■= 
J'^ t and note that in view of the graphical construction specialized for the conditional 
specification (2.6), the relation (2.5) guarantees that the process in some fixed bounded A C R'' 
exhibits exponentially decaying dependencies on the external configuration in D!^ as n ^ 00. 

2.4 Stabilizing functionals of Gibbs point processes 

In this section we specialize to our Gibbs point process setting the notion of a stabilizing functional, 
see [3, 23, 24, 25] and the references therein. As in section 1, let ^(•, •) be a translation invariant 
functional defined on pairs {x,X) where X is a. finite point collection in Mf^ and x G X. Further, 
when a; ^ A", we abbreviate £,{x,X D {x}) by ^{x, X). 

Next, suppose that a given point process S on is stochastically dominated by a homogeneous 
Poisson point process and suppose that there exists C3 > such that for every ball Br{x) the 
conditional probability of Br{x) not being hit by S given the external configuration £ := E\Br{x) 
admits the bound 

¥[E n Br{x) = 0|f ] < exp(-C3r'^) (2.7) 
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unifomily in £. Stochastic domination and (2.7) provide upper and lower stochastic bounds on the 
number of points in any ball analogous to those satisfied by a homogeneous Poisson point process 
and for this reason such S are called Poisson-like. The next proposition tells us that the Gibbs 
point processes considered here are Poisson-like. 

Proposition 2.1 Every Gibbs point process "P* with an exponentially localized potential is a 
Poisson-like process. 

Proof. Indeed, the stochastic domination by V comes from the obvious relation 70 C po in 
the above graphical construction of because po coincides in law with V. The second relation 
(2.7) follows by the graphical construction as well. Indeed, we have A{x, 0) < 00 and hence, by 
(2.2) and (2.3), in the course of the dynamics given by the graphical construction the acceptance 
probability for a birth attempt at some y inside a ball Br-s{x) with no points alive in the whole 
Br{x) is uniformly bounded away from 0, both in the location of the point y attempting to be born 
and in the external configuration, as soon as s and r > s are taken large enough. On the other 
hand, the ball reaches a completely empty state with intensity at most 1. Consequently, the time 
fraction of having the ball fully empty decays exponentially with the volume of the ball uniformly 
in the external configuration and hence so does the probability of having no point alive in Br{x) 
at the time by stationarity of the graphical construction in time. □ 

Similarly to (1.2), say that ^ is a stabilizing functional in the wide sense if for every Poisson-like 
process S there exists an a.s. finite stabilization radius R := R^{x,'E), such that a.s. 

^{x, E n Br{x)) = ^x, [E n Br{x)] U A) (2.8) 

for all locally finite point collections C R'^ \ Br(x). Stabilizing functionals in the wide sense can 
a.s. be extended to the whole process E, that is to say for all a; e M'' 

^{x,E) := lim £,{x,EnBR{x)) 

r—*oc 

is a.s. well defined. 

Given s > and a Poisson-like process E define the tail probability 

sup F[R{x,EnQx}> s], F[R{x,E)>s] . 

Further, we say that ^ is exponentially stabilizing in the wide sense if for every Poisson-like 
process E we have limsup^^,^ log r(s) < 0. Thus, if ^ is exponentially stabilizing in the wide 



t(s) := r(s, E) := max 
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sense, then there exists a C4 such that for all s € M+ we have 

sup P[R{x, E)> s]< exp(-C4s) and sup V[R{x, S n Qa) > s] < exp(-C4s). (2.9) 
xeM<* x>i,xeQx 

We stress that, unlike in the standard Poisson input setting of [3, 23, 24, 25], where Poisson 
points in disjoint sets are independent, the configuration S fi -Bjj(a;) will, in general, depend on 
the configuration in S fi Br{xY. Thus, unlike the standard Poisson input setting, the wide sense 
stabilization of ^ at a; within radius R does not imply that the value of ^(a;, S) does not depend on 
the configuration outside Bii{x); on the other hand this value is independent of the configuration in 
Br{xY given the configuration 'Ef\Bii{x). This weak dependence feature of wide sense stabilization, 
which carries additional technical considerations, allows us to establish limit theory for functional 
and measures in geometric probability over point sets more general than the usual Poisson and 
binomial point sets. 

As we will see shortly, many functionals which stabilize in the standard Poisson input setting 
also stabilize in the wide sense. Possibly there are some functionals which stabilize over Poisson 
samples but which do stabilize in the wide sense, but we are not aware of these functionals. For 
these reasons, when the context is clear, we will henceforth abuse terminology and use the term 
'stabilization' to mean 'stabilization in the wide sense', with a similar meaning for 'exponentially 
stabilizing'. 

2.5 Functionals with bounded perturbations 

The theory presented in this paper is mainly confined to translation invariant geometric functionals 
and its extension to non-translation invariant functionals seems to require non-trivial effort. Nev- 
ertheless, a small step towards the non-translation invariant set-up can be made with only slight 
modification of the existing theory. This extension is the subject of the present subsection and it 
deals with asymptotically negligible bounded perturbations of translation-invariant functionals. To 
put it in formal terms, consider the following notion. Consider a Poisson-like input point process 
S. Assume that £,{■■,■) is a translation invariant geometric functional exponentially stabilizing in 
the wide sense and let ^(•, •; A) be a family of geometric functionals indexed by the extra parameter 
A > 0, not assumed to be translation invariant but enjoying the following properties: 

• For each A > the functional ^(•, •; A) admits a representation 

i{x, X- A) = i{x, X) + 5{x, X- A), (2.10) 
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where the correction (perturbation) 5{x, X; A) is not necessarily translation invariant but, for 
all p > it satisfies the moment bound 

sup E[(5(a;,S;A)]P < £(A,p) < 00, (2.11) 

where lim^^cx) £(A,p) = for each fixed p. 

• The perturbation S{-.-:X) satisfies the wide sense exponential stabilization with the same 
stabilization radius R^{-, •) as ^. 

If these two conditions hold, we say that ^(■, •; A) is an asymptotically negligible bounded perturbation 
of(, on input S,; for brevity we call it just a bounded perturbation of ^ in the sequel. The message of 
this subsection, to be made formal below, is that the asymptotic behavior of bounded perturbations 
of a translation invariant functional is indistinguishable from the asymptotic properties of the 
functional itself. This observation brings the limit theory for stochastic quantization within the 
compass of stabilizing functional ; see section 7. 



3 Weak laws of large numbers and central limit theorems 

We now state our main results, which show that sums of stabilizing functional defined on Gibbsian 
input (with exponentially localized potential) on large cubes satisfy weak laws of large numbers 
and Gaussian limits as the cube size tends to infinity. For all A > 0, let Qx := [-X^/'^/2, X^/<^/2]<^ 
be the volume A cube centered at the origin of M*^, and let //^ be the A-rescaled ^-empirical measure 
on Qi := [-1/2, 1/2]'*, that is 

4- E ?(a^,^*ngA)(5,/Ai/- (3.1) 

xev^nQx 

Let :— /^^(Qi) be the total mass of /x^, and for future reference, define also the non-rescaled 
infinite- volume measure 

/x« := J2 (3.2) 
Let p G [0, oo). Say that ^ satisfies the p- moment condition if 

sup sup £[1^(2;,^* U A')|P] < 00, (3.3) 

where C denotes the collection of all finite point sets in R"*. 
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Recall that B{Qi) denotes the set of bounded / : Qi — > M and that := — ]E[Ma]- 
Under appropriate moment conditions, our first two results establish a weak law of large numbers 
and variance asymptotics for (/, / e B{Qi), as A — > oo. Our third result shows that the 
finite-dimensional distributions of (A~^/^(/i, /U^), A~^/^(/m, /x^)), /i,...,/m G B{Qi), converge 
to those of a multivariate normal as A ^ oo, and, in the univariate CLT we establish a rate of 
convergence. Finally our last general result establishes asymptotics for bounded perturbations of 
a translation invariant ^. 

Theorem 3.1 (WLLN) Assume that S. is stabilizing and satisfies the p-moment condition (3.3) 
for some p> 1. We have for each f G B{Qi) 

lim A-iE[(/,Mi)]-r£(r) / fix)dx (3.4) 

where 

E{t) :=eHt) :=E[e(0,7'*)exp(-A(0,P*))] . (3.5) 
Moreover, if (3.3) is satisfied for some p> 2 then A~^ (/, converges to tE(t) Jq^ f{x)dx in L"^. 

Note that E{t) depends on the underlying intensity t via "P* even though this parameter does 
not explicitly show up in the defining formula. Before stating variance asymptotics write 

a«[0] :=E [e2(0,P*)exp(-A(0,P*))] 

and for all a; G M'' define the two point correlation functions for the functional £, over the Gibbsian 
input 7-** by 

a^[0,x] :=E[e(0,7'*U{x})e(a;,P*U{0})exp(-A({0,ar},O)] -[^^Wr, (3-6) 
where A({a;, y}, X) := *(A' U {x, y}) - *(A'). 

Theorem 3.2 (Variance asymptotics) Assume that ^ is exponentially stabilizing and satisfies the 
p-moment condition (3.3) for some p> 2. We have for each f G B{Q\) 

hm A- iVar[(/, 4)]= ry«(T) / f\x)dx, (3.7) 

where 

V^{t) := a^[0]+T I a^[0,x]dx <oo. (3.8) 
Letting iV(0, tr^) denote a mean zero normal random variable with variance a^, we have: 
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Theorem 3.3 (CLT) Assume that ^ is exponentially stabilizing and satisfies the p-moment con- 
dition (3.3) for some p> 2. We have for each f G B{Qi) 



X-'/'{f,fii)^N{0,rVHr) / f{x)d 



(3.9) 



and the finite- dimensional distributions {X ^/^(/i, /x^), A ^^'^{fm,l^x)), fi, ■-, fm ^ B{Qi), con- 
verge to those of a mean zero Gaussian field with covariance kernel 

(/i,/2)^Ty«(r) / fi{x)f2{x)dx. (3.10) 
JQi 

Moreover, if (3.3) is satisfied for some p > 3 then for all X>2 and all f € B{Qi) we have 



sup 



if, 4) 

Var[(/,/if)] 



< t 



P[7V(0, 1) < t] 



<CilogXf'^X-'/^. 



(3.11) 



Assuming that ^(-j •; A) is a bounded perturbation of a stabilizing functional ^ we have: 

Theorem 3.4 The analogs of Theorems 3.1, 3.2 and 3.3 hold for ^ replaced by £,{-,-;X) under 
their respective stabilization and moment assumptions. 



Remarks, (i) Comparison with [12]. The results of [12] establish limit theory for functionals 
^ of weakly dependent Gibbsian input, but essentially these results require ^ to have finite range 
(finite range test functions). Theorems 3.1-3.4 extend [12] to cases when ^ has infinite range and 
stabilizes. 

(ii) Comparison with functionals on Poisson input. Theorems 3.1-3.4 show that the established 
limit theory for stabilizing functionals on homogeneous Poisson input [3, 21, 23, 24, 25] is insensitive 
to weak Gibbsian modifications of the input. Thus the entirety of weak laws of large numbers and 
central limit theorems for functionals defined on homogeneous Poisson input given previously 
in the literature [3], [21]-[25] extend to the corresponding analogous results for functionals on 
point processes whose local specification (2.1) with respect to the Poisson process is exponentially 
localized. If the input is Poisson, then the term A(a;,P*) vanishes, and hence Theorem 3.1 
extends the Poisson weak law of large numbers in Theorem 2.1 of [25]. Likewise, Theorem 3.2 
extends the variance asympotics of [3] and [21], whereas Theorem 3.3 extends the central limit 
theory of [3], [21] and [26]. 
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(in) Numerical evaluation of limits. We emphasize that the point process is intrinsically 
algorithmic; this algorithmic scheme provides an exact (perfect) sampler [12]. It is computationally 
efficient and yields a numerical evaluation of the limits (3.5) and (3.7). 

(iv) Extensions and generalizations. The variance convergence (3.7) and the asymptotic normal- 
ity (3.9) hold under weaker stabilization assumptions such as power-law stabilization (see Penrose 
[21]), but the resulting additional technical details obscure the the main ideas of our approach, and 
thus we have not tried for the weakest possible stabilization conditions. Similarly, counterparts to 
Theorems 3.1-3.4 should hold for functional defined in terms of non- homogenous Gibbsian input, 
but we do not provide the technical details here either. 

4 Proofs of main results 

This section is organized as follows. First, in Subsection 4.1 we establish exponential clustering 
properties for stabilizing functional of processes with exponentially localized potentials. Expo- 
nential clustering is central to our approach, as it shows that the cumulants of (/, /x^), / e B{Qi), 
converge to those of a normal random variable, that is to say they vanish asymptotically upon 
suitable re-scaling for all orders above two. Then, in Subsections 4.2 and 4.3 we establish Theo- 
rems 3.1, 3.2 and Theorem 3.3 respectively, using either the cumulant techniques developed in [3] 
or the Stein techniques of [26]. Subsection 4.4 provides the proof of Theorem 3.4. 

4.1 Exponential clustering lemma 

Let be a point process with exponentially localized potential and assume that ^ is an expo- 
nentially stabilizing functional in the wide sense. For each Poisson-like configuration S we denote 
by ^[S] this point configuration marked with the values of ^, that is to say each a; e S carries the 
mark ^(a;, S). We have then 

Lemma 4.1 For each k > 1 there exist M > and c := c{k) > such thai for any deterministic 
points Xi, . . . ,Xk G the total variation distance between ^[P*] restricted to the union Bi{xi) U 
... U Bi{xk) and the product of respective restrictions of ^[P^] to Bi{xi), . . . , Bi{xk) does not 
exceed M/c exp(— cmiuj^j dist{xi,Xj)). 

Proof. The statement of the lemma is a consequence of the graphical construction of the process 
■p* and of the exponential stabilization of ^. To see it, observe that the considered total variation 
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distance does not exceed the probability of the event that the random sets 



Ai:=Ano[ (J Br{x)\ 



xe7'*nsi(xi) 



are not all disjoint, where R := R[x,V'^] and where the causal ancestor cone Ario[A] is defined in 
section 2.2. Indeed, if all Ai'a arc disjoint then the values of ^ over all points in balls B\{xi) depend 
on disjoint and hence independent portions of the free birth-and-death process in the graphical 
construction. To complete the proof it suffices now to show that the probability P[Ai n Aj ^ 0] 
decays exponentially with the distance between Xi and Xj for each i and j. Now, this follows because 

• The number of points in n B\{xi) and fl B\{xj) admits super-exponentially decaying 
tails in view of the Poisson domination property of the Poisson-like process P*, 

• For each such point x the stabilization radius R[x, "P*] admits exponentially decaying tails 
by the wide sense exponential stabilization (2.9), 

• Consequently, the diameter of the union UxeP*nSi(xi) ^r{^) ^^ch balls has exponentially 
decaying tails too, 

• Finally, using the exponential decay relation (2.5) for causal ancestor clan sizes in the graph- 
ical construction, we conclude that the diameter of Ai also has exponentially decaying diam- 
eter. 

The proof is hence complete. □ 



There are several ways to prove limit theorems for stabilizing translation invariant functionals. 
To illustrate the new features arising in the setting of functionals of Gibbsian input, we will first 
assume that / is a.e. continuous, that ^ satisfies the moment condition (3.3) for p = 4, and appeal 
to cumulant methods. In this setting we may directly apply the cumulant methods developed in 
Section 4 of [3] (especially those methods used for proving statements (i) and (ii) of Theorem 2.1 
there) and hence we only provide crucial points, referring the reader to [3] for further details. The 
arguments in Section 4 there show that 



4.2 Proof of Theorems 3.1 and 3.2 




(4.12) 
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and 

lim A-^Var[(/,/i|)] = [g(0) + / c{0,x)dx] [ f{u)du, (4.13) 

where the correlation functions c{x), q{x) and c{x, y), x,y & M'', are the respective Radon-Nikodym 
derivatives given by 

E/i^(dx) — c{x)dx, 

E[{ii^{dx)f] = q{x)dx and 

Cov[/i^(rfa;),/i^(dy)] = c{x,y)dxdy, x^y, 

where ji^ is the infinite-volume empirical measure defined in (3.2). Indeed, the main idea as 
briefiy sketched below is to use stabilization, under the guise of the exponential clustering Lemma 
4.1 here, to show that when proving our results, in the A ^ oo limit we can safely replace 
(modulo a correction of order o(A^/^)) the considered expression (/, /j,^) by (./, [/ioo]!) whore 
[Atoo]i := X^xsT'^nQA Now, the last expression coincides with ^Q^fxdfjfi, where 

fx{x) := f{\~^/'^x), X € Q\. Consequently, E(/, = (/, E/x^) for large A is well approximated 
by fxdKjjfi = Jq^ f\{x)c{x)dx = Ac(0) Jq^ f{u)du by translation invariance of /x^ and upon 
a variable change u := X~^/'^x. This yields (4.12) for continuous /. To get (4.12) in the gen- 
eral set-up, that is to say when / e B{Qi) and when ^ satisfies the bounded moment condition 
(3.3) for some p > 1, one can follow the approach of [21]. Likewise, under exponential stabiliza- 
tion, Var[(/, /i^)] is well approximated by Var[/g_^ f\d^^, which by Campbell's theorem, equals 
■^QxxQx ^y)' where Var/x^ := E[/Lt^ — [E/x^] ig) [E/x^] is the variance measure 

of iJ.^; see Section 4 in [3] and references therein for more details on moment measures. Using the 
usual decomposition of the variance measure into the diagonal and off-diagonal component [3] , we 
see that the last expression equals 

/ fl{x)q{x)dx+ f fx{x)fx{y)c{x,y)dxdy. (4.14) 

Using the continuity of f\, the translation invariance of c(x, y) and q{x), and the exponential decay 
of c{x,y) in the distance between x and y as guaranteed by the exponential clustering Lemma 4.1 
with fc = 2, we come to (4.13) as required. To show (4.13) when / G B{Qi) and when ^ satisfies 
the moment condition (3.3) for some p > 2, we may modify the approach of [21]. 

Now, to calculate the correlation fimctions c(-),c(-,-) and q{-) in (4.12) and (4.13), note first 
that, given in M.'^\dx, the probability of observing an extra point of at x is t cxp(— A(x, 'P^))dx 
as determined by the construction of the process in Subsection 2.2, where rdx corresponds to the 
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birth attempt intensity at x whereas exp(— A(a;, "P*)) comes from the acceptance probabihty. Con- 
sequently, Eij.^{dx) = tE^{x,V^) exp{—A{x,'P^))dx and hence 

c{x) = tE^{x, V*) exp{-A{x, 7?*)). (4.15) 

Likewise, 

q{x) = rE,e^(a;,7'*)exp(-A(a;,P*)). (4.16) 

Further, for x,y £ R'', given "P* in K.'^\{dxUdy), the probabihty of observing extra points of P* at 
X and y respectively is cxp(— A({a;, y}, 'P^))dxdy, where again rdx and rdy are the probabilities 
that the birth attempts at x and y were made whereas exp(— A({x, y}, P*)) is the probability that 
they were both accepted. Consequently, 

c{x, y) = r^E^ix, P* U {y})^{y, U {x}) exp(-A({a;, y}, P*)) - c{x)c{y). (4.17) 

In other words, c{x,y) — t^(T^[0, y — x] with cr^[-,-] as in (3.6). The required relations (3.4) 
and (3.7) follow now by putting (4.12) and (4.13) together with (4.15), (4.16) and (4.17) and 
comparing with (3.5) and (3.7). The exponential clustering Lemma 4.1 when combined with the 
moment conditions imposed on ^ implies that the two-point correlation c(.x, y) exhibits exponential 
decay in the distance between x and y whence the integral J^^^ c(0, x)dx is finite. This observation 
allows us to conclude the finitcness of E{t) and V{t), as given by (3.5) and (3.8), respectively. 
Consequently, the L^-convergence stated in Theorem 3.1 follows now by the variance convergence 
in Theorem 3.2 and, given (4.12) and (4.13), the proof of both of these theorems is complete. □ 

4.3 Proof of Theorem 3.3 

When / is continuous on Qi and when ^ satisfies the moment condition (3.3) for all p, the expo- 
nential clustering Lemma 4.1 allows us to use the techniques developed in Section 5 of [3], where it 
replaces the clustering Lemma 5.2, to show that all cumulants of (/, are all of the volume order 
A and hence, upon the A~*^/^-re-scaling with k being the order of the cumulant, the cumulants 
of order higher than two vanish asymptotically yielding the required Gaussian limit; see [3] for 
details. 

More generally, for / G B{Qi) and when ^ satisfies the moment condition (3.3) for all p > 2, the 
rate (3.11) holds by following verbatim the the Stein approach of [26], using wide sense stabilization 
and the exponential clustering Lemma 4.1 instead of stabilization. Combining (3.7) and (3.11) 
yields (3.9) for / e B{Qi). This completes the proof of Theorem 3.3. □ 
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4.4 Proof of Theorem 3.4 

The uniformly decaying bound on all moments of the perturbation term 6{-,-;X) in (2.10), as 
stated in (2.11), combined with the exponential stabilization of the perturbation, allows us to use 
the Holder and Minkowski inequalities to conclude that the addition of 5(-, •; A) to ^ does not affect 
the asymptotic behavior of the first and second order correlation functions. This means that the 
cumulant-based argument for Theorems 3.1 and 3.2 carries over also for ^( ^ •; A) with no further 
modifications. This yields the bounded perturbed versions of Theorems 3.1 and 3.2 for continuous 
test functions and under the moment condition (3.3) with p = 4. To relax the moment conditions 
as in the respective statements and to get the results for general bounded test functions we resort 
again to the approach of [21], which completes the proof of these two theorems for functionals 
with bounded perturbations. The remaining CLT Theorem 3.3 for functionals with bounded 
perturbations follows now easily by the stabilization property imposed on the perturbation term 
S{-, •; A) in full analogy with the respective proof of Theorem 3.3. □ 

5 Examples of Gibbs point processes with exponentially 
localized potentials 

The notion of an exponentially localized potential 'J is general and includes the following non- 
exhaustive list of the corresponding point processes T'*. If an energy functional 4* has finite in- 
teraction range so that its add-one potential satisfies A{x,X) = A{x,X n Br{x)) for some r, as 
would be the case in many examples considered by [12], then clearly (2.4) is satisfied and there 
usually are natural ways of ensuring that the constant C\ is large enough so that exponential 
localization and (2.5) hold as well. These include decreasing the intensity r of the underlying 
Poisson process V which corresponds to increasing A(-, ■) by a positive constant (low reference in- 
tensity/density regime) as well as multiplying ^ and hence also A(-, •) by some small enough /3 > 
(high temperature regime). The following list is not limited to finite range energy functionals. 

(i) Strauss processes. A Strauss process involves perturbing a Poisson process according to an 
exponential of the number of pairs of points closer than a fixed cutoff. For such processes the 
add-one potential depends only on points within the cut-off range and so ip{r) vanishes when r 
exceeds this cut-off. 

(ii) Point processes with pair potential function. A large class of Gibbs point processes [30], 
known as pairwise interaction point processes and including the Strauss process, has Hamiltonian 
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'^{X) := ^ X^i<j — Xj\\), X := {xi}, with (j) bounded below, usually assumed to be positive 
by absorbing the offending constant into the intensity of the underlying Poisson process. If the 
pair potential function ^ has finite range, as would be the case with the Strauss process, then 
the potential ^ is localized since tp vanishes beyond the interaction range. On the other hand, 
suppose the pair potential function has infinite range, but satisfies the following strengthened 
superstability condition: ^ decays exponentially fast and ^{s) = +00 for s < ro, that is there is a 
hard-core exclusion condition forbidding the presence of two points within distance less than ro, 
[28]. In this context then the point process is easily verified to be exponentially localized as 
soon as the intensity r is low enough (low density regime) or (j) admits a sufficiently small upper 
bound on its oscillations (high temperature regime). 

(iii) Area interaction point processes. This is a germ grain process, where the grain shape 
is a fixed compact convex set and where the potential at each Poisson germ is determined by a 
function of the intersection of the grains at that germ. As a special and simple instance, suppose 
that the energy functional £)(<¥) is a scalar multiple 7 of the volume of the union of the radius 
r balls centered at points x £ X (1 D. Then, for 7 small enough, the resulting area interaction 
process (consisting of 'ordered' points for negative interaction parameter 7 and 'clustered' points 
for positive interaction parameter 7) is exponentially localized. More general energy functional 
involve an additive term representing a scalar multiple of the total number of points [1], which can 
be alternatively absorbed into the intensity. As noted in [1], area interaction processes plausibly 
model certain biological processes, including those where the realization of the process represents 
spatial locations of plants (or animals) consuming food within distance r. The energy functional 
is then a scalar multiple of the area of the food supplying region. These are described more fully 
on p. 9 of [12] and in [1]. 

(iv) Point processes defined by the continuum Widom-Rowlinson model. Another example of 
the point process is that defined in terms of the continuum Widom-Rowlinson model from 
statistical physics, see [31] as well as [13]. Here we have fixed radii (say radius equal to a) spheres 
of two types, say A and B, with interpenetrating spheres of similar types but hard-core exclusion 
between the two types. This defines a point process whose potential is exponentially localized as 
soon as the reference intensity is low enough, since the function ^(r) vanishes when r > 2a. It 
is known, see ibidem, that the continuum Widom-Rowlinson admits an equivalent reformulation 
in terms of single-species gas of interpenetrating spheres which is area-interacting in the sense 
of point (iii) above - this is seen by integrating out the positions of B particles and keeping 
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track of the locations of A-particles only. Likewise, upon forgetting the marks carried by the 
particles in the two-species representation one gets the so-called random cluster representation for 
the Widom-Rowlinson model, see [6] and [13], from which the law of the Widom-Rowlinson model 
can be recovered by assigning independent and equiprobable A and B-labels to maximal connected 
clusters of particles, whence the name random cluster model. Theorems 3.1-3.4 are valid for all of 
these equivalent models, provided the intensity is low enough, as discussed above. 

(v) Point processes given by hard-core model. An important and natural model falling into 
the framework of our theory is the so-called hard-core model with low enough reference intensity. 
In its basic version the hard-core model, extensively studied in statistical mechanics, arises by 
conditioning a Poisson point process on containing no two points within distance less than 2r 
for some r > standing for a parameter of the model. Clearly, this process admits a Gibbsian 
description with ^ set to +oo if there are two points closer than 2r from each other and otherwise. 
Consequently, the potential is exponentially localized if the reference intensity of the underlying 
Poisson point process is low enough or if r is small enough, which also reduces to decreasing the 
reference intensity upon appropriate re-scaling (in fact rather than imposing separate conditions on 
r and the reference intensity r it is enough to require that rr'^ be small enough, as easily checked 
by re-scaling). 

(vi) Truncated Poisson process. The hard-core gas is a particular example of a truncated Poisson 
process. In general, a truncated Poisson process arises by conditioning a Poisson point process on 
the event that a certain family of constraints is fulfilled. In this paper the constraints imposed 
are of the following form: we fix a certain family of bounded sets and require that none of these 
sets contain more than a certain given number of points. Such processes are used in modelling of 
communication networks [2]. In particular, if we require that no ball of radius r contain more than 
some constant number k of Poisson points, then ip vanishes beyond r and the associated point 
process has an exponentially localized potential, possibly upon decreasing the intensity. 

6 Applications 

Below we indicate some applications of our main results. This list is not exhaustive and does 
not include applications to e.g. germ-grain models where the germs arise as the realization of the 
Gibbsian point process with an exponentially localized potential. 
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6.1 RSA packing with Gibbsian input 

Let X c 'R'^ he locally finite. Consider a sequence of unit volume d-dimensional Euclidean balls 
Bi,B2,---- with centers arriving sequentially at points in X. The first ball Bi to arrive is packed and 
recursively, for i = 2,3, ... let the ith ball be packed if it does not overlap any ball in Bi,B2, ....Bi_i 
which has already been packed. Let ^{x, X) be either or 1, depending on whether the ball arriving 
at X is either packed or discarded. 

When X is the realization of a Poisson point process on Q\, this packing process is known 
as random sequential adsorption (RSA) with Poisson input on Q\. When X is the realization of 
an infinite sequence of independent random d- vectors uniformly distributed on the cube Q\, then 
this is called the RSA process with infinite binomial input; in such cases, RSA packing terminates 
when it is no longer possible to pack additional balls. In dimension d = 1, this process is known 
as the Renyi car parking problem [27]. In the infinite input setting and when d = 1 Renyi [27] 
(respectively Dvoretzky and Robbins [9]) proved that the total number of parked cars satisfies a 
weak law of large numbers (respectively central limit theorem) as A — > oo; recently these results 
were shown to hold for all dimensions in [21] and [29]. 

Virtually all limit results for RSA packing assume that the input is either Poisson or a fixed 
number of independent identically distributed random variables. To the best of our knowledge, 
RSA packing problems with Gibbsian input have not been considered before in the literature. The 
following theorem widens the scope of the existing limit results for RSA packing. Put 

A~ E ^(^'^*nQ;,)J,/;,i/., 

SO that N{'P'^ n Qx) := J2xev'^nQx ^(^' ^ '3^) denotes the total number of balls packed on Qx 
from the collection of balls with centers in "P* n Qa- 

Theorem 6.1 Let be Gibbsian input with an exponentially localized potential. Then 

A-i7V(7?* n Qx) ^ tE [e(0, P*) exp(-A(0, P*))] in (6.1) 

and 

A-iVar[Ar(P* n Qx)] ^ tV^{t) (6.2) 

where V^{t) is given by (3.8). The finite- dimensional distributions (A~^/^(/i, /x^), \~^^^{fm, mI))) 
/ij •••) fm € B{Q\), converge to those of a mean zero Gaussian field with covariance kernel 

{fi,f2)^TVHr) [ h{x)f2{x)dx. (6.3) 
JQi 
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Remark. As spelled out in [24], Theorem 6.1 also applies to related packing models, including 
spatial birth growth models with Gibbsian input as well as RSA models with balls replaced by 
particles of random size/shape/charge, and ballistic deposition models. 

Proof. The approach used in [24] shows that the packing functional ^{x, •) is exponentially 
stabilizing on Poisson-like sets. Indeed, any Poisson-like set S can be coupled on the common 
underlying probability space with a dominating Poisson point process Vt of finite intensity r, r 
large, and containing S a.s. Now, the idea underlying the argument in [24] shows that the packing 
status of a point a; in a configuration X depends on X only through its algorithmically determined 
sub-configuration CZ[a;, Af] referred to as the causal cone or causal cluster of x in the presence 
of X, see [24] for details. The causal cluster Cl[x,X] is easily seen to be non-decreasing in X. 
In particular, using that S C yields CZ[a;,!E!] C Cl[x,Vr] a.s. for x € E. However, by the 
arguments in section 4 of [24], the causal clusters generated by points of Vr exhibit exponential 
decay, and hence so do causal clusters of points in S showing that the packing functional ^{x,-) 
is exponentially stabilizing on Poisson-like sets, in particular on P*. In other words, ^(x, ■) is 
exponentially stabilizing in the wide sense. Clearly ^ satisfies the bounded moment condition (3.3) 
and therefore Theorems 3.1- 3.3 show that the {f,fi\), f S B{Qi), satisfy the weak law of large 
numbers and central limit theorem given by (6.1- 6.3), respectively. □ 

6.2 Functionals of Euclidean graphs on Gibbsian input 

In many cases, showing exponential stabilization of functionals of geometric graphs over Poisson 
point sets [3, 23], can be reduced to upper bounding the probability that regions in R'^ are devoid of 
points by a term which decays exponentially with the volume of the region. When the underlying 
point set is Poisson, as in [3, 23], then we obtain the desired exponential decay. When the underlying 
point set is Poisson-like, the desired exponential decay is an immediate consequence of condition 
(2.7). In this way the existing stabilization proofs for functionals over Poisson point sets carry 
over to stabilizing functionals on Poisson-like point sets. This extends central limit theorems for 
functionals of Euclidean graphs on Poisson input to the corresponding central limit theorems for 
functionals defined over Gibbsian input. The following applications illustrate this. 

(i) fc-nearest neighbors graph. The fc-nearest neighbors (undirected) graph on the vertex set 
X, denoted NG{X), is defined to be the graph obtained by including {x, y} as an edge whenever y 
is one of the k nearest neighbors of x and/or x is one of the k nearest neighbors of y. The fc-nearest 
neighbors (directed) graph on X, denoted NG'{X), is obtained by placing a directed edge between 
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each point and its fc-nearest neighbors. 

Total edge length of k-nearest neighbors graph. Let L{X) denote the total edge length of NG{X) 
and let ^(a;, X) denote one half the sum of the edge lengths of edges in NG{X) which are incident 
to X. Put 

Theorem 6.2 Let "P* be Gibbsian input with exponentially localized potential. Then 

A-iL(P*nQA)^TE[e(0,7'*)exp(-A(0,7'*))] m L'' (6.4) 

and 

A-iVar[i(P* n Qx)] ^ rV^ir) (6.5) 

where V^{t) is given by (3.8). The finite-dimensional distributions {X~^^'^ {fi, ft^^) X~^/'^ {fm, fix)), 
fi, •••) fm G !3{Qi), converge to those of a mean zero Gaussian field with covariance kernel 

{fl,f2)^TV^{T) [ fl{x)f2{x)dx. (6.6) 

JQi 

Remark. Theorem 6.2 generalizes Theorem 6.1 of [23], which is restricted to nearest neighbor 
graphs defined on Poisson input. 

Proof. Let S be a Poisson-like point set. Considering the arguments in the proofs of Theorem 
6.1 of [23] and Theorem 3.1 of [3], it is easily seen that the set of edges incident to any point 
X in NGCE) is unaffected by the addition or removal of points outside a ball of random radius 
R. Moreover, the radius R has exponentially decaying tails, which may be seen as follows. For 
simplicity we prove exponential stabilization in dimension two, but the argument is easily extended 
to higher dimensions by using cones instead of triangles (for d = 1 we use intervals instead of 
triangles). For each t > construct 6{k + 1) triangles Tj{t), 1 < j < 6{k + 1), such that a; is a 
vertex of each triangle and such that each triangle with edge containing x has length t. Let Rx 
be the minimum t such that each triangle contains at least one point from S. In such a situation, 
the union of the 6(A; + 1) triangles Tj{t), 1 < j < 6{k + 1), may be partitioned into 6 equilateral 
triangles with common edge length t, each triangle containing at least k + 1 points. Then, because 
S is Poisson-like, it follows that F[Rx > t] < 6{k + 1) exp(— Csf"^). Moreover, as explained in the 
proof of Theorem 6.1 of [23], simple geometry shows that AR^ is a radius of stabilization for the 
functional (, at x. Thus ^ is exponentially stabilizing. 
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An easy modification of the proof of Lemma 6.2 of [23] shows that L moreover satisfies the 
p-moments condition (3.3) for all p. Therefore Theorems 3.1- 3.3 show that the (/, f £ B{Qi), 
satisfy the weak law of large numbers and central limit theorem given by (6.4- 6.6), respectively. □ 

Number of components in nearest neighbors graph. Let k = 1. Given a locally finite point 
set A", let Sy\x,X) denote the reciprocal of the cardinality of the component in NG{X) which 
contains x. Thus H{X) := J2xex ^'''^^i^^ denotes the total number of components of NG{X). 
Put 

Theorem 6.3 Let be Gibbsian input with exponentially localized potential. Then 

\-^H{V^ n Qx) ^ rE [^(0, P*) exp(-A(0, P*))] in (6.7) 

and 

A-iVar[if (P* n Qx)] ^ tV^{t) (6.8) 

where V^{t) is given by (3.8). The finite- dimensional distributions (A^"^/^(/i, /l^), {fm^ A*a))' 
/i) ■■■jfm G l3{Qi), converge to those of a mean zero Gaussian field with covariance kernel 

{h,h)^TV^{r) [ fi{x)f2{x)dx. (6.9) 
Jqi 

Proof. We establish that ^1^1 is exponentially stabilizing on Poisson-likc sets S and appeal to 
Theorems 3.1-3.3. When k = 1, the Poisson-likc properties of the input process and the methods 
of Haggstrom and Mecstcr [16] and Kozakova, Mccster. and Nanda [18] show there are no infinite 
clusters in A^G'(S). Moreover, the proof of Theorem 1.1, 1.2 and Proposition 2.2 of [18] and property 
(2.7) of Poisson-like processes show that the finite clusters in NGC^) have (super)exponentially 
decaying cardinalities and diameters. Now we show exponential stabilization of ^ as follows. Let 
Rx be the radius of the cluster in iVG'(S) containing x. Write Br for Br{0). Put 

R := sup R^. 

xeBn^ns 

It is not hard to see that R has exponentially decaying tail. Indeed, writing 

oo 

P[R > f] = V P[ sup R^ > t, 2-^-1 <Ro< 2^] 

OO OO 

< ^ P[ sup i?x > t, 2-^-^ < Rq] < ^exp(-C2^)P[ sup R^ > t] 
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and then noting that the cardinaUty of x G B23 fl S decays polynomially fast in 2^ with overwhelm- 
ing probabiUty, we obtain the desired exponential decay of R. We can also show that AR is a radius 
of stabilization for ^I'^l at the origin (see proof of Lemma 6.1 of [23]) . Since ^I'^l trivially satisfies the 
bounded moments condition (3.3) for all p, the weak law and central limit theorem for H{V^ nQx) 
and A-^/2(/i,/i^) follows by Theorems 3.1- 3.3. □ 

(ii) Voronoi tessellations. Given X cM."^ and x G X, the set of points in closer to x than 
to any other point of A" is a convex polyhedral cell C(x, X). The collection of cells C{x, X), x € X, 
form a partition of R'' which is termed the Voronoi tessellation induced by X. 

Total edge length. Given X C M.'^, let L{x, X) denote one half the total edge length of the finite 
edges in the cell C{x, X). It is easy to see that L is exponentially stabilizing on Poisson-likc sets S. 
Indeed, when d = 2, it suflttccs to follow the arguments in the proof of Theorem 8.1 of [23] and to 
note that stabilization radius depends on finding a minimum edge length t such that 12 isosceles 
triangles with this edge length have at least one point from S in them. Because S is Poisson-likc 
we may follow the arguments in [23] verbatim to see that L stabilizes. See section 6.3 of [21] for 
the case d > 2. 

As in section 6.3 of [21], we may show that L also satisfies the moment condition (3.3) for 
p = 3. It follows that the total edge length X^seP^nQA L{x,V^ C\ Q\) of the Voronoi tessellation 
on HQa satisfies the weak law and central limit theorem as A ^ 00. In other words, the Voronoi 
analog of Theorem 6.2 holds. 

(ill) Other proximity graphs. There are further examples where showing exponential sta- 
bilization of functionals of geometric graphs (in the wide sense) involves upper bounding the 
probability that regions in are devoid of Poisson-like points. Such estimates are available in 
the Poisson setting and it is not difficult to extend them to Poisson-like point sets. 

In this way, by modifying the methods of [23] (sections 7 and 9) and [3] (section 3.1), we obtain 
weak laws of large numbers and central limit theorems for the total edge length of the sphere of 
influence graph, the Delaunay graph, the Gabriel graph, and the relative neighborhood graph over 
Gibbsian input P*. 

6.3 Gibbsian continuum percolation 

Let be a locally finite point set and connect all pairs of points which are at most a unit distance 
apart. The resulting graph is equivalent to the basic model of continuum percolation, in which 



24 



one considers the union of the radius 1 balls centered at points of X, see Section 12.10 in [15]. Let 
^'■^l {x, X) be the reciprocal of the size of the component containing x, so that 

iV(A') = ^e[^](x,A') 
xex 

counts the number of finite components in G. 

Section 9 of [23] discusses central limit theorems for N{T' n Qx). Using Theorem 3.3 we can 
generalize these results to obtain a central limit theorem for the number of components N{V^ f\Qx) 
in the continuum percolation model on Gibbsian input in the subcritical regime, possibly of interest 
in the context of sensor networks on Gibbsian point sets. To formulate this central limit theorem 
assume that our Gibbs point process with exponentially localized potential is such that it 
admits a stochastic upper bound by a homogeneous Poisson point process of some intensity r 
falling into the subcritical regime of the considered continuum percolation (see Section 12.10 in 
[15]). Note that due to the Poisson-like nature of the input process this is always possible 
upon spatial re-scaling. We argue that ^I'^l is exponentially stabilizing on Poisson-like sets S 
as follows. If S is Poisson-like and if r is subcritical, then S is also subcritical by stochastic 
domination. Consequently, the diameter of the connected cluster emanating from a given point 
has exponentially decaying tails, see ibidem. This yields the required exponential stabilization 
upon noting that ^^'^^{x, ■) does not depend on point configurations outside the connected cluster 
at X. Moreover, ^1^1 is bounded above by one and thus satisfies the moments condition (3.3). Hence 
by Theorems 3.1-3.3, NCP^OQx) satisfies the weak law of large numbers and central limit theorem, 
exactly as in the statement of Theorem 6.3. 

6.4 Functionals on Gibbsian loss networks 

Given the Poisson point process V Vr, consider the following Gibbs point process. Fix an 
integer m G N. Attach to each point of 7^ a bounded convex grain K and put the potential to 
be infinite whenever the grain K at one point has non-empty intersection with more than m other 
grains. This condition prohibits overcrowding, and, for more general repulsive models, one can 
put large and finite whenever the grain K at one point has non-empty intersection with a large 
number (some number less than to) of other grains. The resulting point process, which we call "P*, 
represents a version of spatial loss networks appearing in mobile and wireless communications. As 
discussed in point (vi) in Subsection 5 the so defined \1/ is exponentially localized as soon as the 
underlying intensity r is small enough. 
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Let K be an open convex cone in M'' (a cone is a set that is invariant under dilations) with 
apex at the origin. Given x,y G , we say that y is connected to x, written a; — > y, if there is 
a sequence of points {xi}"^^ G {fC + x) (1 V^, \xi — Xi+i\ < 1, \xi — a;| < 1 and \y — Xn+i\ < 1- 
If the length of this sequence does not exceed a given m, we write x — V- For all r > let 
B!^{x) ■.= x + {ICnBr{0)). 

Coverage functionals. The functional ^(a;, P*) := sup{r e M : x — > y for all y € B^{x)f\V^} 
determines the maximal coverage range of the network at x in the direction of the cone K. 

The coverage measure is := X^j;£p*nQx ^(^'^* ^ Qa)^x/Ai/<*- Confining attention to V, 
where r belongs to the the subcritical regime for continuum percolation, is in turn subcritical 
because of Poisson domination. Since the continuum percolation clusters generated by any Poisson- 
like set S have exponentially decaying diameter, it follows that ^ stabilizes in the wide sense 
(recall the proof for the number of components in the continuum percolation model) and that ^ 
admits an exponential moment. By appealing to Theorems 3.1 and 3.3, we obtain a weak law of 
large numbers and central limit theorem for both the coverage measure and the total coverage 
ExeP-nQ.^(^>^*nQ,). 

Network reach functional. Say that the network has reach at least r at a; if a; — > y for all 
y G B^{x) n P*. Put ^riXjV^) := 1 if the network has reach at least r at a; and otherwise put 
^r{x, ■P*) := 0. Theorems 3.1 and 3.3 yield a weak law of large numbers and central limit theorem 
for the total network reach J^xev'^nQx ^r{x,V^ fl Q\). 

Number of customers obtaining coverage. Independently mark each point x of "P* with mark T 
(transmitter) with probability p > and with mark R (receiver) with the complement probability. 
Then define the reception functional ^(a;,^*) to be 1 if a; is marked with T or (when x is marked 
with R) a z ^ X for some z in the transmitter set {z e : z marked with T}. Put ^(a;, P*) 
to be zero otherwise. Thus ^(a;, •) counts when a customer at x gets coverage and the limit theory 
for the sum X^a:e'P*nQA ^* ^ Qx)^ which counts the total number of receivers (customers) 
obtaining network coverage, is given by Theorems 3.1 and 3.3. 

Connectivity functional. Given a broadcast range r > and the transmitter set {z G : 
z marked with T}, let ^r{x,V^) be the minimum number, say m, such that every point in y e 
B^{x)(lV^ can be reached from some transmitter z G with m or fewer edges or hops, that is to 
say there exists a transmitter z such that z -^m V for all y G B^(x) nP*. Thus all receivers in the 
broadcast range r > can be linked to a transmitter in m or fewer hops. Small values of £,r{x, V^) 
represent high network connectivity; for each r > 0, Theorems 3.1 and 3.3 provide a weak law 
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of large numbers and central limit theorem for the connectivity functional '^^r{z,V^ <^Q\) as 
A — > oo. 

7 Gibbsian quantization for non-singular probability mea- 
sures 

Quantization for probability measures concerns the best approximation of a d-dimensional prob- 
ability measure P by a discrete measure supported by a set Xn having n atoms. It involves a 
partitioning problem of the underlying space and it arises in a variety of scientific fields, including 
information theory, cluster analysis, stochastic processes, and mathematical models in economics 
[14]. The goal is to optimally represent P, here assumed non-singular with density h, with a point 
set Xn, where optimality involves minimizing the stochastic quantization error (or 'random 
distortion error') given by 



Recall that for all x and locally finite point sets X, C{x,X) denotes the Voronoi cell ('Voronoi 
quantizer') induced by the Euclidean norm around x with respect to X. 

The optimal (non-random) quantization error is given by min;t^ I{Xn) and the seminal work 
of Bucklew and Wise [5] shows that this minimal error satisfies 



where ||/i||d/(d+r) denotes the d/{d + r) norm of the density h and where the so-called rth quanti- 
zation coefficient Qr,d is some positive constant not known to have a closed form expression. 

The first order asymptotics for the distortion error on i.i.d. points sets (that is to say letting Xn 
consist of i.i.d. random variables) was first investigated by Zador [32] and later by Graf and Luschgy 
[14] and Cohort [7]. Letting X„ be i.i.d. random variables with common density h'^/('i+'-) / J 
and Wd the volume of the unit radius d-dimensional ball in R'^, Zador's theorem shows 




lim n''''^m.mI{Xn) = Qr,d\\h\\d/{d+r) 



(7.1) 



lim n''/'^I{Xn) 



co^'^^^Til + r/dmid/^d+r) 



(7.2) 



■oo 



whence (see Prop. 9.3 in [14]) the upper bound 



Qr,d < l^d 



-r/d- 



T(l+r/d). 



(7.3) 
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Molchanov and Tontchev [20] have pointed out the desirabihty for quantization via Poisson 
point sets and our purpose here is to estabhsh asymptotics of the quantization error on Gibbsian 
input. This is done as follows. For A > and a finite point configuration X we abbreviate 
X(^X) '■= A~^/'^Af. Moreover, we write X := X f) Qi so that in particular := X~^/'^X n Qi. 
Recall also that we write for as in the previous sections. Consider the random point 
measures induced by the distortion arising from V^^-^ , namely 



Wc will be interested in the asymptotic behavior of the random integrals (/, /i*). Clearly, 
when / = 1 then (/, /j*) gives another expression for the distortion /(P*^). On the other hand, 
if / = 1b, then (/, ■) mcasm-es the local distortion. This section establishes mean and variance 
asymptotics for (/, /i*} as well as convergence of the finite-dimensional distributions of (/, /i*). 
Since we will no longer be working with translation invariant ^, we will need to appeal to Theorem 
3.4. 

Put 

M*(r) / |M;rdwexp(-A(0,7'*)) 

where, recall, t is the intensity of the reference process V. Note that M*(t) does depend on r 
through 7-"*. Changing the order of integration we have 

EM*(r)=E / |y ripens,,, (y)=0exp(-A(O,P*))rfj/= (7.4) 



/ 



|y rE [exp(- A(0, P*))lp*ns,„ iy))=0]dy. 

In the special case where = (i.e. coincides with the reference process V) and where 
the intensity r of P is 1 we readily get EM°(1) = r(l + 3)0;^''^'*. More generally EM°(r) = 

^-(l+r/d)r(l + L)u;-^/'i. Put 

F*(t) :=E[M*(r))2] + 



/ (E / \w\''dw [ \w-y\''dwexp[-A{{0,y},P'^) 

JR<* V JC(0,P*U{j/}) Jc{v,V^U{0}) 



(EM*(r))2 Uy. 



For any random point measure p, recall that p denotes its centered version, that is p := p — Ep. 
Theorem 7.1 Assume that the density h of P is continuous on Qi. We have for each f £ B{Qi) 
Urn X'-/'^{f,lil} =tE[M'^{t)] [ h{x)f{x)dx inL^ (7.5) 
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and 

lim Ai+2''/'^Var[(/,/iJ)] =rF*(T) / f{x)h'^{x)dx. (7.6) 

The finite- dimensional distributions X^^^^^^^'^({fi,Jlf), (/fc,/!*)), /i, ...,/fc G of the ran- 

dom measures (A^^/^+'"/''7I*) converge as X ^ oo to those of a mean zero Gaussian field with 
covariance kernel 

(/i,/2) ^ Ty*(r) / fiix)f2ix)h\x)dx, fuf2 e (7.7) 

Upper bounds for quantization coefficients. When / = 1 the right hand side of (7.5) 
gives 

Hm y/'^{l^nl)= hm A''/''/(:P(*,) =tE[M*(t)]. 

Since the Bucklcw and Wise limit (7.1) is necessarily no larger than the right hand side of the 
above, this shows that in addition to the bound (7.3), that the rth quantization coefficient Qr,d 
also satisfies the upper bound 

Qr,d< (||/l|U/(d+r))-VE[M*(T)]. 

Recall from our discussion above that when ^E* = (i.e. is Poisson) and when / = 1, then the 
right hand side of (7.5) equals t~'^I'^u:'2^^'^V{\ + r/d) and thus 

Qr4 < m\d/(d+r))-'T-^/v^''n-i-+r/d). 

We believe, although are not yet able to provide a full proof, that whereas the distortion error 
(7.5) is relatively large for Poisson input, it can be made smaller if we restrict to point sets which 
themselves enjoy some built-in repulsivity while keeping the same mean point density. Indeed, 
given a fixed mean number of test points it seems more economical to spread them equidistantly 
over the domain of target distribution than to allow for local overfuUs of test points in some regions, 
which only result in wasting test resources with the quantization quality improvement considerably 
inferior to that which would be achieved should we shift the extraneous points to regions of lower 
test point concentration. In other words, the right hand side of (7.5) for repulsive Gibbs point 
processes should be smaller than the corresponding distortion for the Poisson point process with 
the same point density. It should be emphasized here that in order to stay within the set-up of our 
asymptotic theory we have to assume that the repulsivity is weak. On the other hand, it is very 
likely that going to some extent beyond this requirement may lead to even smaller quantization 
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errors. These seem to be natural and interesting questions, yet at present we cannot handle them 
with our current techniques. 

Proof of Theorem 7.1. We claim that the assertions of Theorem 7.1 can be reduced to an 
application of Theorem 3.4 for functional with bounded perturbations. We do it first assuming 
that the density h is bounded away from 0. To this end, consider the following parametric family 
of geometric functional : 

,,P(dA-i/<^y) 



JC{x,X) 



Xh{X-^/dx) 



\y-x\ \\-i/d \ dy. 
c(x,x) h{\ Vda;) 



Putting 



i{x,X):= f \y-x\^dy 

JC(x,X) 

we obtain the bounded perturbed representation (2.10) for ^(•, •; A) with 



6{x,X;\) :-- 



C(x,X) h{X l/da;) 



(7.8) 



(7.9) 



(7.10) 



It is easily seen that on Poisson-like input both ^ and d as given in (7.9) and (7.10) stabilize 
exponentially with common stabilization radius determined by the diameter of the Voronoi cell 
around the input point; see Subsection 6.2(ii). We claim that this 6{-, •; A) also satisfies the bounded 
moments condition (2.11). To this end, use the fact that h is bounded away from and write 



\d{x,X;X)\ < C 



f 

JC{x,X) 



y-xY{h{X-"-'y)-h{X-"''x))dy. 



(7.11) 



Letting X := and using the exponential decay of the Voronoi cell diameter on Poisson-like 
input, see Subsection 6.2, we conclude from (7.11) that, for all }3 > and A > 

p 



supE[5(a;,P*;A)f < supCE 



/ \y-x\-{h{X-^/''y)-h{X-^/''x))dy 

Jcix.V^) 



lC{x,p-^) 

By the translation invariance of and by the uniform continuity of the density h, this is bounded 
above by 



CE 



C(o,-p*) 



\yruJhiX-'/%\)dy 



=■■ L{p,r, A), 



is dominated by 



where ujh{') is the modulus of continuity of h. Since /^j^q p^,^ \y\^aJh{X ^'"^{yDdy 
an intcgrablc function of lu uniformly over A, namely by a constant multiple of the p{r + dY'^ power 



of the Voronoi cell diameter, and since 



/c(o,P*) \yrM>^-'/''\y\)dy 



converges to zero as A — > oo 



for almost all ui, we may use the bounded convergence theorem to conclude that L{p,r,X) 
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as A ^ 00. This is clearly enough to get (2.11) and hence ^(-j ■; A) is an asymptotically negligible 
bounded perturbation of ^(-j •). To proceed, we note that the quantization empirical measure /x* 
satisfies for each / e B(Qi) 

{f,^^t) = X-'-'/''{fh,4), (7.12) 
where is the standard empirical measure (3.2) for | as in (7.9), that is to say 

On the other hand, it is easily verified that ^ satisfies all assumptions of our limit Theorems 3.1, 
3.2 and 3.3. Consequently, Theorem 3.4 can be applied for ^, which yields Theorem 7.1 via the 
formula (7.12) allowing us to translate results for /i| to the corresponding results for /U*. This 
completes the proof of Theorem 7.1 for h bounded away from 0. 

To proceed, assume now that h fails to be bounded away from and, for e > put := 
max(/i,e) and let /i*^ be the version of /i* with h replaced by h^. Using the definition of /i*, 
and the exponential decay of the diameter of Voronoi cells in a Poisson-like environment we easily 
conclude that 

|E[(/,M?>-(/,mU]I = 0(A-'-/'^£). (7.13) 
Likewise, using the same we get 

Var[(/,M^) - = 0{X-'-'^/''e). (7.14) 

Applying Theorem 7.1 for hg, which is legitimate due to he being bounded away from 0, and 
then using (7.13) and (7.14) we readily get the required expectation and variance asymptotics for 
(/, /U*) as well as the weak law of large numbers, which follows by the variance convergence. 
The remaining central limit theorem statement for (/, /i*) follows directly by the Stein method as 
in Theorem 3.3, which is not affected by h being not bounded away from 0. This completes the 
proof of Theorem 7.1 for general h. □ 
Acknowledgements. J. E. Yukich thanks Ilya Molchanov for helpful discussion related to quan- 
tization. 

References 

[1] A. J. Baddeley and M. N. M. van Lieshout (1995), Area interaction point processes. Ann. Inst. 
Statist. Math., 47, 4, 601-619. 



31 



[2] F. Baccelli, B. Blaszczyszyn, and M. Karray (2006), A spatial markov queuing process and its 
applications to wireless loss systems, preprint. 

[3] Yu. Baryshnikov and J. E. Yukich (2005), Gaussian limits for random measures in geometric 
probability. Ann. Appl. Prob., 15, lA, 213-253. 

[4] Yu. Baryshnikov, P. Eichelsbacher, T. Schreiber, and J. E. Yukich (2008), Moderate devia- 
tions for some point measures in geometric probability. Annales de I 'Institute Henri Poincare - 
Probabilites et Statistiques, to appear. 

[5] J. A. Bucklew and G. Wise (1982), Multidimensional asymptotic quantization theory with rth 
power distortion measures. IEEE Trans. Inf. Th. 28, 239-247. 

[6] J. T. Chayes, L. Chayes and R. Kotecky (1995), The analysis of the Widom-Rowlinson model 
by stochastic geometric methods. Comm. Math. Phys. 172, 551-569. 

[7] P. Cohort (2004), Limit theorems for random normalized distortion. Ann. Appl. Prob., 14, 
118-143. 

[8] D. J. Daley and D. Vere- Jones (1988), An Introduction to the Theory of Point Processes. 

Springer- Verlag . 

[9] A. Dvoretzky and H. Robbins (1964), On the "parking" problem. MTA Mat Kut. Int. KzL, 
(Publications of the Math. Res. Inst, of the Hungarian Academy of Sciences), 9, 209-225. 

[10] R. Fernandez, P. Ferrari and N. Garcia (1998), Measures on contour, polymer or animal 
models. A probabilistic approach, Markov Processes and Related Fields 4, 479-497. 

[11] R. Fernandez, P. Ferrari and N. Garcia (2001), Loss network representation of Ising contours. 
Ann. Prob. 29, 902-937. 

[12] P. Ferrari, R. Fernandez and N. Garcia (2002), Perfect simulation for interacting point pro- 
cesses, loss networks and Ising models. Stock. Proc. Appl. 102, 63-88. 

[13] H.-O. Georgii, O. Haggstrom and Ch. Maes (2000), The random geometry of equilibrium 
phases. In C. Domb, J. Lebowitz (eds) Phase transitions and critical phenomena 18, 1-142, 
Academic Press, London. 

[14] S. Graf and H. Luschgy (2000), Foundations of Quantization for Probability Distributions. 
Lecture Notes in Mathemmatics, vol. 1730. 



32 



[15] G. Grimmett (1999), Percolation (2nd edition). Grundlehren der Mathematischen Wis- 
senschaften, 321, Springer. 

[16] O. Haggstrom and R. Mccstcr (1996), Nearest neighbor and hard sphere models in continuum 
percolation. Random Structures and Algorithms 9, 295- 315. 

[17] P. Hall (1988), Introduction to the Theory of Coverage Processes. Wiley. 

[18] 1. Kozakova, R. Meester, and S. Nanda (2006), The size of components in continuum nearest- 
neighbor graphs. Ann. Prob. 34, 528-538. 

[19] I. Molchanov (1997), Statistics of the Boolean Model for Practitioners and Mathematicians. 
Wiley. 

[20] I. Molchanov and N. Tontchev (2007), Optimal Poisson quantisation, Stat. Probab. Letters 
77, 1123-1132. 

[21] M. D. Penrose (2007), Gaussian hmits for random geometric measures, Electronic J. Prob., 
to appear. 

[22] M. D. Penrose (2007), Laws of large numbers in stochastic geometry with statistical applica- 
tions. Bernoulli, 13, 4, 1124-1150. 

[23] M. D. Penrose and J. E. Yukich (2001), Central limit theorems for some graphs in computa- 
tional geometry. Ann. Appl. Prob. 11, 1005-1041. 

[24] M. D. Penrose and J. E. Yukich (2002), Limit theory for random sequential packing and 
deposition. Ann. Appl. Prob. 12, 272-301. 

[25] M. D. Penrose and J. E. Yukich (2004), Weak laws of large numbers in geometric probability. 
Ann. Appl. Prob. 13, 277-303. 

[26] M. D. Penrose and J. E. Yukich (2005), Normal approximation in geometric probability. In 
Stein's Method and Applications, Lecture Note Series, Institute for Mathematical Sciences, 
National University of Singapore, 5, A. D. Barbour and Louis H. Y. Chen, Eds., 37-58, Elec- 
tronically available via http : //www . lehigh . edu/ ~ j eyO/publications . html. 

[27] A. Renyi (1958), On a one-dimensional random space-filling problem. MTA Mat Kut. Int. 
Kzl., (Publications of the Math. Res. Inst, of the Hungarian Academy of Sciences) 3, 109-127. 



33 



[28] D. Ruelle (1970) Superstable interactions in classical statistical mechanics. Comm. Math. 
Phys. 18, 127-159. 

[29] T. Schrcibcr, M. D. Penrose, and .J. E. Yukicli (2007), Gaussian limits for multidimensional 
random sequential packing at saturation. Comm. Math. Phys., 272, 167-183. 

[30] D. Stoyan, W. Kendall, and J. Mecke (1995), Stochastic Geometry and Its Applications. John 
Wiley and Sons, Second Ed. 

[31] B. Widom and J.S. Rowlinson (1970), New model for the study of liquid- vapour phase tran- 
sitions. J. Chem. Phys. 52, 1670-1684. 

[32] P. L. Zador (1982), Asymptotic quantization error of continuous signals and the quantization 
dimension. IEEE Trans. Inf. Th., 28, 139-149. 

Tomasz Schreiber, Faculty of Mathematics and Gomputer Science, Nicholas Gopernicus Uni- 
versity, Torun, Poland: tomeks@mat.uni.torun.pl 

J. E. Yukich, Department of Mathematics, Lehigh University, Bethlehem PA 18015: 
j oseph.yukich@lehigh.edu 



34 



